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ABSTRACT 

A highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova SN2006aj has 
puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object 
were put forward to explain its multiwavelength signature. We propose that the source is naturally explained 
by another channel, a tidal disruption of a white dwarf (WD) by an intermediate mass black hole (IMBH). The 
tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some 
debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. 
We successfully fit soft X-ray spectrum with the Comptonized blackbody emission from a jet photosphere. The 
optical/UV emission is consistent with self-absorbed synchrotron from the expanding jet front. The accretion 
rate temporal dependence M{t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH 
mass is found to be about IQ/^Mq in three independent estimates: (1) fitting tidal disruption M{t) to soft X-ray 
lightcurve; (2) computing the jet base radius in a jet photospheric emission model; (3) inferring the central BH 
mass based on a host dwarf galaxy stellar mass. The supernova position is consistent with the center of the host 
galaxy, while low supernova ejecta mass is consistent with a WD mass. High expected rate of tidal disruptions 
in dwarf galaxies is consistent with one source observed by Swift satellite over several years at GRB060218 
distance of 150 Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth. 
Subject headings: accretion - black hole physics - gamma rays: bursts - radiation mechanisms: general - 
supernovae: general - X-rays: individual (GRB060218) 



1. EVITRODUCTION 

The existence of the stellar mass black holes (BHs) with 
mass Mbh < IOQMq and the supermassive black holes 
(SMBH) with mass Mbh > lO^M© has long been established. 
A population of intermediate mass black holes (IMBH) likely 
exists. Those BHs have masses in a range IQQMq < Mbh < 
IO^Mq. They could hve in the centers of dwarf galaxies 
( |Dong et al.|[200 7', 'Greene"2012') or globulai- clusters (GCs) 
iFabbiano et al.| 1 997, 2001 ; Colbert & Mushotzky 1999; Mat- 
Isumoto et alpOOlj |Gultekin et al.||2004|). The I MBHs can 
form by a collapse of a massive cloud dBegelman et al. 2006) 
or a massive star ([Fryer et al . 2001; Madau & Rees 2001] 
Schneider et al.|200^ i or grow from a stellar mass BH. A star 
cluster Gl in M3 1 gala xy is estimated to host a2x IO'^MqBH 



(Gebhardt et al. 2602) based on the velocity dispersion pro- 
lile. SDSS J160531.84H-174826.1 dwarf galaxy is estimated 
to have an IMBH with a mass ^ 7 x IO^M q based on lumi- 
nosity scaling relations (|Dong et"aL 2007]l. One of the best 



43-49 galaxy. Based on 



candidates is HLX-1 source in ESO : 
a thin disk thermal state the estimated BH mass is 3 x lO^M© 
( [Davis e t al.||2011|l. All those candidates are still tentative 
and the uncertainty of the mass estimates is up to an order of 
magnitude. More IMBH candidates with qualitatively differ- 
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ent observational signatures may provide stronger evidence 
for existence of such objects. 

Tidal disruptions of the white dwarfs (WDs) by IMBHs pro- 
vide such qualitatively different signature. Unlike disruptions 
of main sequence (MS) stars, which happen far away from 
the BHs and have slow timescales, the disruptions of WDs 
are very fast. Su ch disruptions ma y lead to accretion rates 
up to lO"^M0yr"^ ( Haas et al.| 2012 1. A flow with an extreme 
accretion rate may produce ashort powerful burst of radia- 
tion, for example, when a jet is launched. The X-ray emission 
of Swift J 1644-1-57 source is attributed to a jet launched by a 
super-Eddington accretion disk formed after a tidal disruption 
([Bloom et al.|[20Tll [Burrows et al.[|20TT] [Levan et al.[[2011| 
[Zauderer et al.|201 l| i. 

When particle acceleration is inefficient, then the only 
radiati on from the jet is the photospheric emis s ion ([Eich- 



, ler & Levinson[|2000| [Meszaros & Rees| 120001 [Daigne & 
|Mochkovitch[[2(J02t [Rees & Meszaros[[2005| l. The radiation 
field is in a thermal equilibrium with matter within a dense 
hot jet up to the distance Rp^ along the jet, called a photo- 
spheric radius, where the optical depth to Compton scatter- 
ing is about unity t„ ^ I. The photospheric emission has 



a quasi-blackbody spectrum (Pe'er et al. 2007 Pe'er 2008 
Pe'er & Rydel[20n] [Beloborodov[[201 Ij . Part cle accelera- 
tion may happen either in the entire volume of the jet or in 
the internal shocks. Dissipation in the volume of the jet, in 
particular, in the sub-photospheric regions leads to a modi- 
fied blackbody spectrum from near the je t photosphere ( [Pe'erj 
[eraL][2005l 0OE[ [Gianniosl[2006l |20n| . Modificationsli> 
elude the Comptonization and the broadening of the black- 
body peak with pair produ ction. Effective particle accelera- 
tion in the internal shocks ([Meszaros|2006|) is likely respon- 
sible for the emission of gamma-ray bursts (GRBs). When 
two shells moving at different speeds collide above the pho- 
tosphere, a shock forms, and a substantial fraction of rela- 
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live kinetic energy can be transferred int o the electrons (Rees 
& Meszaros||1994[ [Sari & Piran||1997| l, which radiate syn- 



chrotron emission. ITie collisions of shells are especially effi- 
cient, when the jet bulk Lorentz factor F is large. For the small 
Lorentz factors and the large photospheric radii, such as in the 
jets from tidal disruptions of the WDs by the IMBHs, the col- 
lisions of shells may happen within the photosphere, and the 
internal shock signatures might be weaker The fluctuations 
of the jet F diffuse out, while the jet travels within the pho- 
tosphere, and the fluctuations outside of the photosphere are 
small. Even when the slow shells collide, the energy release is 
weak. Therefore, we expect to see the strong blackbody sig- 
nature of the photospheric emission. While the energy density 
of the slow tidal disruption jet might not warrant fast pair pro- 
duction, the Comptonization is still expected to modify the 
spectrum. 

A distinct feature of the WD disruptions by IMBHs is a su- 
pernova. As a consequence of the tidal compression along the 
angular momentum axis , the WD may undergo thermonuclear 
ignition (Lumin etl 19851 |Luminet & Marck|1985[ |Luminet & 
[Kchon 1989b ) and explode. The energy and the composition 
of the resulting supernova are the functions of the W D mass 
niwD^ the pericenter radius Rp, and the BH mass M^h ( |Ross- 
Iwog et al. 2009). If the WD is massive and the disruption is 
deep, then the explosion could be similar to a supernova Type 
la with a comparable energy release. Less massive WDs may 
lead to low luminosity explosions with little Ni synthesized. 
There is no explosion at all in some cases. In any case the 
ejected mass < IMq is less than in the core-collapse super- 
novae. 

In sum, we predict a transient similar to a GRB, but softer 
and longer, accompanied by a supernova with a small ejecta 
mass. The best candidate we find in Swift GRB catalogue is 
GRB060218 source. This unusual event is an underluminous 
very long GRB with a duration fgo ~ 2600 s and with a smooth 
X-ray lightcurve (Soderberg et al. 2006). It is accompanied 
by a fast supernova SN2006aj, whic h was modeled to hav e 
a low ejecta mass Mej ^ (1 -2)M0 (MazzaU et al. 2006b| ). 
The X-ray emission has a blackbody component characteris- 
tic of the photospheric radiation. The early X-ray radiation 
is accompanied by the powerful optical/UV emission (Ghis- 
[ellini et al. ,2007a,) . Two classes of theories were propo sed to 
explain the source: a supernova shock breakout model (Cam- 



pana et al.|2006}|Waxman et al.[2007]jNakar & Sari|2012| ) and 



a model with a mildly relati vistic jet from a newborn compact 
object, such as a magnetar (Fa n et al.|2006[[Tbma et al.|2007[ 
IGhisellini et al. [200 7b a). Nevertheless, we think that a tidal 
disruption of a WD by an IMBH is not only a viable model 
for the source, but also more naturally explains some features, 
such as the duration and the soft quasi-thermal spectrum. 

The paper is organized as follows. In § [2] we describe in 
more detail the observations of GRB060218/SN2006aj. We 
extensively discuss in §|3]the former theoretical explanations 
of the source: the shock breakout model and the jet launched 
by the newborn compact object. Then in §[4] we concentrate 
on modeling within a tidal disruption scenario. We rederive 
the dynamics of a tidal disruption and discuss the physics 
of jet launching. We perform time-resolved spectroscopy of 
the source and successfully model the X-ray emission with a 
Comptonized blackbody spectrum from the photosphere. The 
derived jet base radius corresponds to about lO^M© IMBH. 
We propose that the origin of the powerful early optical/UV 
radiation is the front region of propagating jet. We find that 



the full X-ray lightcurve can be fitted well by a scaled de- 
pendence of the fallback accretion rate on time M{t). The fit 
provides the IMBH mass about We find that the steep 

decay phase is consistent with the action of absorption alone 
and that the afterglow can be naturally explained as powered 
by the central engine operating in a shallow regime. The 
associated supernova is consistent with a WD origin, while its 
position is consistent with the center of a host galaxy. The 
mass of the dwarf host galaxy provides the estimate of a cen- 
tral IMBH mass about lO^M©. In §[5] we estimate the event 
rates. A high rate of tidal disruptions is predicted in the dwarf 
galaxies, while the disruptions of the WDs constitute a signifi- 
cant percentage of all disruptions. Since the disruptions of the 
MS stars mostly happen at large distances from the BHs, then 
such events could be much dimmer, than the disruptions of 
the WDs. The tidal disruptions of the WDs provide a plenty 
of material to feed the central IMBH up to a supermassive 



2. OBSERVATIONS OF GRB060218 
GRB060218 triggered the Burst Alert Telescope (BAT) on- 



board Swift mission satellite on 18 Feb 2006 (Campana et al. 
[2006). Soon after the trigger the X-ray Telescope (XRT) 
mentified a bright source, whose count rate peaked at around 
960 s from BAT initial trigger. Then the source gradually de- 
cayed over a continuous observation period, which ended at 
around 2700 s. The event duration tgo = 2600 s is unusually 
long among Swift GRBs (Campana et al. 2006). The XRT 
lightcurve was unusually smooth and regular. The rise period 
is characterized by hard emission with the BAT flux about 
equal to the XRT flux (i Campana eFall [20061 |Toma et al.] 
2007 1 and with the peak energy Ep ^ 5 keV. The isotropic- 
equivalent luminosity is Liso ^ 10"^^ er g s~\ which i s about 10^ 
times less than for the typical GRBs ( Toma et al.|2007) . The 
source is found to have a blackbody component with a low 
temperature T ~ 0.2 keV, whose flux contribution increased 



with time from 15% tiU 80% ( [Campana et al. 2006; Li 2007 ). 
The Swift satellite returned to observe the source at ^ 6000 s 
and found it in a steep decay phase with the absorbed flux 
~ 100 times below the peak value. After about 1 .5 x 10"^ s the 
decay of the X-ray flux flattened into a shallow afterglow with 
the X-ray luminosity Lx (x (Soderberg et al. 2006). The 
simultaneous observations by the Ultra-Violet/Optical Tele- 
scope (UVOT) onboard Swift satellite revealed substantial 
early emission peaking at about 5 x 10"^ s with a deredde ned 
UV flux i^F,, - 10" "'erg s'^cm'^ ( [Ghisellini et al.|2007a| . A 



lower peak at a level I'Fi, « 6 x 10 erg s 
ter several days. The second peak is attributed to a supernova 
SN2006aj CMaeda et al."2007''Ma zzaU et al.|2 006b!l Modjaz 



et al. 2006; Pian et al. 2006 ), while the origin o f the first peak 
is debated ( Waxman et al.|2007i|Ghisellini et al.|2007b]al). 

The supernova was classified as Type Ic ( [Mazzah et aL] 
2006b ). However, its unique spectral propert ies prompted to 
suggest a new Type Id classification (Mazzali et al. [2006b 



Maeda et al. 2007). The rapid supernova peaked at 10 days 
(jPianet al. 2006 ), which is the fastest of all supemo vae asso- 
ciated with GRBs ( [Ferrero et al.|2006| ). According to |MazzaIi] 
( 2007 ), "the eje cted mass pre dicted below 8000km s~ 



et al. 



by the model used by [Mazzali etal., ( ^2006b^ is = IM^ 



while the nebular phase observations suggest a higher total 
ejecta mass Mej = 2Mq ( [Mazzah et"aL[[2006b[ [2007[ [Maeda 



et al.|2007 ). The estimated ^''Ni mass in the ejecta is O.2M0. 
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The ejecta consists mostly of oxygen and carbon (|Mazzali| 
et al.||2007[l, but also contains some silicon and iron (Mira- 
The correspondent supernova energy is 
Ek ^ 2 X 10^'erg. Both the ejecta mass and and total en- 
ergy are much less than the typical values Me, = IOMq and 



Eg = 3 X 10^'erg for GRB-supemovae (see 



Mazzali et al. 



2006b| and references therei n). Substantial optical polariza- 
tion (I Gorosabel et al.||2006"l l indicates the asymmetry of the 
ej ec ted material. 
The optical observation s of the host galaxy indicate a sourc e 



redshift z = 0.0335 ( Mirabal et al.|2006JjMod]az et all 

which corresponds to a distance d = 143 Mpc. GRB060218 
is associated with a dwarf star- forming ga laxy with the stellar 
mass Msi ^U^-^Mq (Ferrero et al. 2007 \ the metalHcity Z ss 
Q.QIZq ( jWierse ma et al. 2007 ), and the characteristic radius 
/?8o = 0.55 kpc (Svensso n et al. 2010) . The supernova and 
the host galaxy were observed by Hubble Advanced Camera 
for Surveys (ACS) instrument under the program GO 10551 
(PI Kulkarn i) in cycle 14. The images produced by [Misra 
et al. ( 201 1[ ) reveal some irregularity of the host dwarf galaxy 
morphology. 

3. FORMER THEORETICAL MODELS 

GRB060218 and the associated supernova SN2006aj 
sparked substantial interest among the researchers with over 
ten papers being dedicated to the theoretical explanations of 
this highly unusual source. The theoretical efforts can be di- 
vided into two big categories: a shock breakout model and a 
model with a jet launched by a magnetar or by a stellar mass 
BH. In the following subsections we briefly review the theo- 
retical models emphasizing their successes and problems. In 
the following section we offer an explanation of the source 
within the tidal disruption scenario. 

3.1. Shock Breakout Model 

The shock breakout model states that the nuclear explo- 
sion following the collapse of a massive core launches a 
shock wave, which propagates out through the star (Colgate, 
[1968 ). The radiation-dominated shock dep osits its energy mto 



me low-density gas at the stellar surface ( |Katz et aL 2010| l. 
The heated gas radiates bremsstrahlung photons, which are 



g PH 

inverse-Compton scattered into the X-ray and the 7-ray bands. 
Gas exponentially cools down. The shock deposits up to 
lO^^^e rg of energy, most of which is radiated as an X-ray 
flash ([ Colgate|1974| l. A variety of sources were successfully 
modeled with the supernova shock breakout. A good exam- 
ple is XRF080109 and the associated Type Ib/c supernova 
SN2008D. The large supernova ejecta mass for that source 
Mej = 4-7Mo (Soderb erg et aL)|2008| [Mazzali et al.||2008|i 
and the low X-ray radiation energy 

10fc _i04bgj.g ( |chevalier| 
|& Fransson|2008| l leave little doubts about its shock breakout 
origin. Other source s include the X- ray brightening of SNLS- 
04D2dc supernova (Schawinski et al. 2008) with low X-ray 
radiated energy < 10^'erg, SN1998bw with the ejecta mass 



12M q accompanied by XRF/GRB980425 ( |Woosley et aL 
1999]), SN20031W with the ejecta mas s I3M0 acc ompanied 



5>rrweak XRF/GRB031203 ( MazzaUet al.|2006a| ). 



In the case of XRF060218 the deposited energy lO'^^'^erg 
is quite large (Campana et al. 2006). The peak isotropic un- 
absorbed flux of the soft X-ray component is Fbb = (3 - 6) x 
10"^ erg s~'cm ~^ with the temperature Tbb = 0.11-0.17 keV 
(|Butler||2007|). These parameters correspond to blackbody 



(3-8) X 10'^ cm, while [Campana et al.| ( |2006| inferred R = 
(0.5-1) X 10'^ cm within their analysis at the early times and 
R = (0.2-2) X 10'^ cm at the late times. This radius indi- 
cates a compact progenitor such as a Wolf-Rayet (WR) star. 
Large explosion energy might be inconsistent with the non- 
relativistic expansion. Mildly relativistic eject a velocity with 
V » 0. 85c was proposed in later modeling by |Waxman et al. 
( |2007[ ), who also stated that the emission should come from 
the photosphere above the stellar surface. The relativistic 
shock breakout theory was further developed by Nakar & SariJ 
( [2012) , who computed the shock dynamics, the lightcm^ve ana 
the spectrum. The spectrum is pred icted to be quasi-the rmal 



with the temperature around 50 keV. |Nakar & Sari[ ( 2012]) de- 
rived a relation between the XRF event duration, the total en- 
ergy, and the observed temperature 



tto 20s 



/ Eho 



\ lO^^erg 



1/2 



'ho 



50keV 



-2.68 



(1) 



To satisfy this relation the early temperature of 40 keV is 
taken for GRB060218 consistent with the peak energy 36 keV 

at the very early times ( |Toma et al.|20 07i 

The shock breakout theory has its problems. |Li|(j2007)) con- 
cluded that the observed temperature and the total energy of 
GRB060218 lead to the unrealistically large photospheric ra- 
dius Rph of a WR star inconsistent with the galactic WR stel- 
lar population. This conclusion is strengthened by our larger 
estimate of the emission radius. The large R^^ in the model 
by Nakar & Sari (2012 ) corresponds to unrealistically dense 
stellar wind. Gh isellini et al., (2007 b) showed that the opti- 
cal/UV spectrum of the source is too bright for the X-rays and 
the optical to be the parts of the sa me blackbody component. 
In response, Waxman et al. (2007 > "fiercely argued'(^that the 
shock breakout is anisotropic, and the different regions emit at 
the different wavelengths, which helps to reconcile the theory 
with the observations. The prompt thermal X-rays are emitted 
by a compressed shell, while the optical radiation originates 
in the outer shells of the e xpanding star at a much larger ra- 
dius. In their final reply, Ghisellini et al.|( |2007a[ ) found the 
anisotropic expansion unconvincing. The simple energetics 
argument is dramati cally inconsisten t with the shock break- 
out model by jWaxman et al. ( 2007 ). The isotropic optical 



thermal emission at t ~ 10^ s requires the temperature about 
^ 1 X 10'' K, which corresponds to the total energy carried by 
a shock above 10^' erg, much in excess of any shock breakout 
model. 

3.2. Jet Launched by a Magnetar or a Stellar Mass BH 

A promising alternative explanation for GRB060218 is the 
jet launched by the central engine. A newborn central engine, 
a magnetar or a BH, results from a core collapse. The jet 
pierces through the star and escapes. That is how a typical 
GRB operates and our source can be just on a low-luminosity 
end of the spectrum. A jet with a low Lorentz factor F ^ 5 and 
a wide opening angle 0^0.3 can produce the emi ssion, which 
peaks in the X-rays and exhibits no jet break (Toma et al. 



2007 1. Such a jet may be powered over a long timescale by a 



magnetar. The central engine activity is the natural explana- 
tion of the afterglow (Soderberg et al. 2006; Fan et al |2006| ). 
[Dainotti et al.| ( [2007) proposed a model of the interaction of a 
electron-positron fireshell and circumburst medium to explain 



emission from the non-relativistic medium at a radius R = 



Ghisellini et al. 2007a 
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the afterglow. Their best model constrains the density profile 
to be « cx r"^, where /3 = 1.0-1.7 up to 10"^ cm. The un- 
usually bright prompt optical/UV emission can be explained 
as produce d by the self-absorbed Comptonized synchrotron 
( [GhiseUini et al.||2007b[ l, while the blackbody X-rays resul t 
from the photospher ic jet emission ( Ghisellini et al.|[2007a| l. 
[Bufano et al.| ( |2012 l argued that the central magnetar would 
spin down too rapidly and could not power the source over 
thousands of seconds, while expending mo st of its energy in 
the exp ansion of the dense envelope. Yet, jQuataert & Kasen 
(|2012!) proposed a model, in which a magnetar can power a 
GRB source over many days or even weeks. 

The relativistic jet is a natural explanation for the source. 
However, if GRB060218 was a GRB, then it clearly was 
an unusual one. The long GRBs have large total energies 
Eiot iso = 10^^"^^ erg (Nava et al. 2008 ), short durations fgo = 
2 - 200 s (Butler et al. 2007 ), an d~fiigher Lorentz factors 
r = 100-1 000 ( ,Pe'er et al.„2007i |Liang et al.|[20T0l [Ghls^ 
ellini|2012) l. The high Lorentz factor of a typical long GRB 
is the consequence of jet confinem e nt by the pressure from 
the star (Tchekhovskoy et al.||2010p. [ Bromberg et al.|(|2011)l 
argued that the GRB population does not extend to the slow 
low-luminosity end. Thus, another origin of events such as 
GRB060218 is to be sought. Below we propose that the jet 
is launched by a tidal disruption. This model provides the ex- 
planations for the observed event duration, the spectrum, and 
the lightcurve. 

4. TIDAL DISRUPTION MODEL 

GRB060218 can be best modeled by a low-luminosity wide 
mildly relativistic outflow/jet. The tidal disruption of a WD 
by an IMBH produces such a jet. In this section we discuss 
the expected temporal and spectral properties of such tidal dis- 
ruptions, fit observations of GRB060218, and consider the ac- 
companying supernova SN2006aj and the host galaxy. Exten- 
sive observations of the source pose multiple tests for any the- 
ory aimed to explain it. The tidal disruption scenario passes 
all these tests. 

4.1. Disruption Dynamics 
4.1.1. Fallback Material and Disk Formation 

A Star gets tidally disrupted, when it approaches sufficiently 
close to the BH. Let us define a tidal radius as 



Rt 



f2M, 



BH 



1/3 



(2) 



V '"WD 

Only the stars on orbits with the pericenter dis tance Rp 



smaller than about Rp < Rj can get disrupted (see Evans & 
[Kochanek ,1989. a nd references therein). The WD radius is 
( |NauenbergfT972 1 



-1/3 



/ mwD \ 



4/3- 



1/2 



(3) 



where Mch = 1.44M0 is the Chandrasekhar mass. We cali- 
brated the normalization to reproduce a typical observed WD 
radius (Nalez yty & Madej 2004) for one solar mass niwo = 
IMq. a tidal radius depends weakly on the BH mass, while 
the event horizon radius Rg = GMbh/c^ is proportional to the 
BH mass. Thus, the SMBHs will swallow the WDs with- 
out disruptions and onl y the BHs with re latively small masses 
can disrupt the WDs ( [Luminet & Pichon ^1989a^ ^Rosswog, 



et al. 20091. For a star on a parabolic orbit about 50% of 
Its material ends up being bound to the BH and about 50% 
remains unbound in a standard picture (Evans & Kochanek 
|1989| l. The fractions may change, when the pericenter dis- 
tance is close to the BH marginally bound orbit radius (several 
gravitational radii). Among mo dels with different pericenter 
radii, BH spins, and orientations Haas et al. ( 20121 found the 
cases, when almost all material is captured or almost all ma- 
terial remains unbound. 

The bound material mak es its way onto the BH. In a classic 
theory the fallback time is ( |Evans & Kochanek|1989| l 



Rl 



VCM^Ri' 



3/2- 



(4) 



The fallback time tf\ , can be any where from several hours 
down to 3 min (Haa s et aL]|2012 1 for the disruptions of the 
WDs by the IMBHs. The peak accretion rate is achieved at 



^peak ~ ^ 



Rl 



Vgm^r 



3/2 ' 



(5) 



where S ~ 3.33 for the adiabatic i ndex Fad = 5/3 based on 



simulations by ,Evans & Kochanekj ( ,1989) and Laguna et al. 
(1993) for Pt =Rt /Rp = F while^'^5.5 for Fad = F4 
( jLodato et al.|2009] l. The smaller values of the adiabatic in- 
dex Fad < 5/3 are more appropriate for the heavy WDs. The 
correspondent peak fallback rate is 



■Mpeak ~ O.OSniwD 



Rl 



,3/2 



(6) 



for Fad = 1 -4. The fallback rate rises from zero, when the most 
bound debris just reach the BH, up to the maximum value 
Mpeak at fpeak and then decreases according to t~^^^ law. 

The fallback matter forms an accretion disk at a circulariza- 
tion radius R^rc = several x Rp. The disk exists as a radiatively 
inefficient accretion flow (RIAF), since th e material density is 
very high and the photons cannot escape (Abramowicz et al. 
|1988| l. The fallback rate determines the accretion rate during 
the early evolution of the system. After the disk is formed, 
it starts to evolve on a slow viscous timescale for Rfirc radius 



( Cannizzo et al.(T990| . Such evolution changes the te mporal 
slope of the accre tion rate to f~^/^ for adiabatic RJAFs ( [Can- 
nizzo et al.|201 l| l. RIAFs with the out ward energy flux may 



unbind the ma terial in the outer disk ( Narayan et al.|[2001[ 
|Metz ger|2012|l. Then only a small percentage of the material 
reaches the BH. However, the inner disk regions with the ra- 
dius less than ^ 2QRg evolve fa st and are not influenced by 
relatively slow energy transport ( Abramowicz et al.|2002[ l. 

As was recently shown by general relativistic magneto 
hydrogynamic (GRMHD) s imulations ( [Tchekhovskoy et al.| 
[201 1 1 McKinney et al.[2()T2| i, a poloidal magnetic field is nec- 
essary toTiuncE~ajerTlie field strength for a strong jet can 
be estimated based on the equipartition argument. When the 
magnetic field energy density is comparable to the internal 
energy density of plasma as 



(7) 



where Wnuc is the density of nucleons, nip is the proton mass. 
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Then the magnetic field near the event horizon is 



B 



BH 



1.3 X 10 



10 / 



M 



1/2 



( 



BH 



V lO^Mgyr-i 
The magnetic field can reach 10'^ 



VlO^Mg 



G. (8) 



a GRB jet (|Usov 
|MacFady"en| 



2007 



10 "' G at the base of 

|T992| [Wheeler et al.||2000[ [Uzdensky &| 



ITaMwaki et al.|[2009| l.^ powerful jet 

is launched after a WD disruption by an IMBH for Bgff ^ 
10'" G. How can such a strong poloidal magnetic field be pro- 
duced? 

4.1.2. Generation of Magnetic Field 

The magnetic field lines get advected towards the BH event 
horizon with the fallback material . The magnetic field of a 



WD is typically within B < lO"* G (|Angel 1978 Putney 



1999: 



with only a few examples of a stronger field. Since the WD 
radius is on the order of the BH gravitational radius, then no 
field amplification is expected due to compression in the con- 
verging fallback flow. The BH may possess an accretion disk 
with a large magnetic field even before the disruption hap- 
pens. When the accretion rate is close to the Eddingto n rate, 
the BH sustains the so-called Eddington magnetic field (Rees 
p84, Daly 2011. ) 



BEdd = 6x10' 



/ Mb 



-1/2 



G. 



(9) 



The magnetic field on the order of B ^ 6 x 10'' G may already 
exist near the horizon of a lO'^M© BH. This field is higher 
than the advected magnetic field of a WD, but is still lower 
than the equipartition magnetic field for the peak accretion 
rate following a tidal disruption. 

In situ amplification of the initial seed magnetic field into 
the equipartition poloidal magnetic field is needed to sustain a 
powerful jet. The magneto-rotational instability (MRI) (Bal- 



bus & Hawley 1991 1998) can exponentially increase the 



magnetic field strength. The amplification rate of the turbu- 
lent magnetic field (to be distinguished from the instability 
growth rate) corr esponds to a local e-folding time of 3 orbital 
periods ( [Stone e t al. 1996) under the MRI action. When start- 
ing from the Eddington field strength, the magnetic field needs 
to grow 10^ times to reach the equipartition. The growth by a 
10^ factor can be achieved over ~ 25 local orbits. The growth 
timescale is short, when the most bound material settles near 
the innermost st able circular o rbit (ISCO) of t he BH. The 
orbital period is ([Bardeen et ar.[|1972i [Shapiro & Teukolsky[ 

\m6j — 



tnrh — 



2ttR^ 



3/2 



+ a 



2.5 



/ Mb 



\ 10* Mr: 



(10) 



at the ISCO for the dimensionless spin a* = 0.6. In sum, a 
IQ'^Mq BH can amplify the random magnetic field 10^ times 
over as little as 60 s. 

The MRI turbulence produces the random B-field, while the 
regular poloidal magnetic field is needed to launch a jet. The 
magnetic field generati on in a turbulent medium s hould occur 
via the dynamo action ( [Brandenburg & Subrama nian 2005 ), 
which operates on a viscous timescale determined by t he ra- 
dial velocity v, ( [Davis et al.|2010| [O'Neill et al.|2011| . The 
geometrically thick accretion flow with the eff ective dimen- 
sionless viscosity a = 0. 1 -0.3 ( (King et al.|2007| l has the radial 



velocity Vr < av^, where vk is the Keplerian velocity. Then 
the viscous timescale at the ISCO can be as short as 15 s for 
a 1O'*M0 BH. Substantial generation of the regular poloidal 
field component is expected over > 10 viscous timescales or 
> 150 s. The failure of the previous attempts to generate the 
magnetic field via the dynamo action might arise from the low 
resolution in those simulations (Jonathan McKinney, private 
communication). In sum, it might be possible to generate the 
substantial poloidal magnetic field before the estimated accre- 
tion rate peak time. 

4.1.3. Jet Launcliing 

The jet is a relativistic outflow of the material from near 
the BH. The spinning BH surrounded by the accreting mag- 
netize d gas is expected to launch a jet (Blandford & Znajek] 
[1977| l. We can estimate the jet power via a Blandford-Znajek 
formula. When the magnetic field energy density reaches the 
equipartition with the matter energy density close to the BH 



then the Blandford-Znajek jet power is approximately ( McK- 
inney|2005||Tchekhovskoy et al.|2010| 



'Poa*^Mc^. 



(11) 



While the simulations with the maximum magnet i c flux were 
able to reac h Po ^ 1 ( Tchekhovskoy et aLl|2011| [McKinney 
[et al.[pOT2] l, more traditional values observed in the simula- 
tions with t he weak in itial field are Pq = 0.01-0.1 ( McKinney] 
& Blandford 2009; P enna et al.|2010| ). The tidal disruption of 
a WD by an IMBH has a typical maximum accretion rate of 
^peak ^ lO'^M^yv'^ according to the formula (j6|. This cor- 
responds to the jet power Lkin ^ O.OlMc^ = 5 x lO'^^^erg s"'. 
Since there is no theoretical understanding whether or not the 
dynamo action would quickly generate the poloidal magnetic 
field, the kinetic power estimated by the equation (fTTl) is an 
upper limit. The actual jet power may be substantiallylower. 

The kinematics of the jet is characterized by the bulk 
Lorentz factor F. The GRB jets have very high F = 100- 
1000 likely because of con finement by the surrounding star 
( [Tchekhovskoy et al. [[2010 1). In turn, the active galactic nu- 
clei (AGN) jets are not strongly confi ned by the ambient 
gas. Their F-fac tors are about F ^ 10 (Jorstad et al. 2005 
Pushkarev et al.[2009). Similarly, the jets from the low-mass 



X-ray b inaries (LM>^Bs) have b ulk Lorentz factors about 
F - 10 ( [MiUer-Jones et al [[2006| . The tidal disruption de- 
bris is expected to scatter in radius, cool, and provide little 
pressure support for the jet regardless of the direction, where 
the debris is scattered. Therefore, the tidal disruption jet is ex- 
pected to have F ^ 10. A low Lorentz factor F ^ 10 is indeed 
suggested for Swift J1644-H 57 tidal disruption event ( |MetzgCT| 
[etal.|2012[[Liuetal.|2012| ). 

4.2. Prompt X-ray Emission 

Prompt X-ray emission from a slow jet is expected to 
have a Comptonized blackbody spectrum produced near the 
jet photosphere. We perform time-resolved spectroscopy of 
Swift XRT observations of the source and fit the spectrum 
with th e Comptonized blackbody model. We follow Butler] 
( [2007[ ) and cut the full XRT window mode observations into 
11 time slices with about 16,000 photons in each. We run 
xrtpipeline, select the time slices with xselect, take the ap- 
propriate response matrix files (RMF) from the calibration 
database (CALDB) (version 1 1) as indicated by xrtmkarf rou- 
tine within the pipeline, and use the ancillary response files 
(ARE) generated by xrtpipeline. We run grppha to group 
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nearby bins to have at least 20 photons per bin and per- 
form the minimization of least squar es. We model the spec- 
trum in XSPEC vl2.7 ( [ArnaudligW ) with the broken power- 
law bknpower and with the blackbody spectrum Compto nized 
by thermal electrons as a part of a compPS model (Pouta- 



nen & Svensson J996|^ . We favor the exact compPS model 



over the approximate prescriptions for thermal Comptoniza- 



tion o f the blackbody spect rum such a s compBB (Nishimura 
et al.| ri986) and compTT ( |Titarchuk| [T994). We compute 



the model flux in (0.05 -10) ke V band to capture ve ry soft 
thermal X-rays. Unlike l Butler| ( [2007 ) and Campana et al.| 
(12006), we consider the host gala xy to be a low -metalhcity 



absorber with Z = 0.07 Zq (jWiersema et al.|2007 1, while fol- 



lowing [Campana et al. ( 2005 1 we hx the CJalactic column at 
Nh = 9 Ay. lO'^^cm""^ with the solar metallicity. We ignore 
any potential changes of the hydrogen absorption column with 
time. We search for a joint best fit to all time slices with a sin- 
gle host galaxy Nh- 

The results of spectral fitting are presented in Table [T] for 
the broken power-law fits. A substantial degeneracy exists 
between the soft power-law slope in a bknpower model and 
the host Nh- We fix the soft photon index to be Fi = -1, 
which is representative of the Rayleigh-Jeans tail /v « 
of the blackbody. The host galaxy hydrogen column Nh = 
0.750 X lO^^cm"^ provides the best joint fit to all time slices 
with = 5463.0 for dof = 4836 degrees of freedom. The 
break energy is 0.7-0.8 keV and depends weakly on the to- 
tal flux. Since the emission becomes substantially softer with 
time, the absorbed flux constitutes a much smaller fraction 
of the unabsorbed flux at the late times. This masks the true 
temporal evolution of the source flux. 

The results of spectral fitting are presented in Table [2[ for 
the blackbody spectrum Comptonized by thermal electrons 
as computed within a compPS XSPEC model. The optical 
depth of Compton scattering is fixed at = 1 for scatter- 
ing at the photosphere. The host galaxy hydrogen column 
Nh = 1.090 X 10^'cm"^ provides the best joint fit to afl time 
slices with = 5358.9 for dof = 4836 de grees of freedom. 
The fit is slightly better for this model compared to the bkn- 
power model. Our fit has a higher blackbody flux frac tion at 
the peak (40%) c ompared to ^ 15% fraction in a fit by [Cam- [ 
[pana et ar[ (12006 ), who employed the solar metallicity of the 
absorber and used a sum of the power-law and the blackbody. 
Table [2[ also sum marizes application of the fireball model 
( Pe'er et al.|2007) to the blackbody emission component. We 
can estimate the bulk Lorentz factor of the flow F and the 
jet base radius knowing the temperature Tbb, the observed 
isotropic flux Fbb, the distance to the source, and assuming 
a fiducial ratio Y = 10 of the total fireball energy to the energy 
emitted in X-rays/7-rays. The obtained jet parameters are 
surprisingly consistent between different time slices, which 
might indicate that the fireball model captures the physics 

of the event. The bulk Lorentz factor of F « 2.1Y^q^ corre- 
sponds to a mildly relativistic jet, while the jet base radius is 

Ro = (1.0-1.7) X IO^'^Y'q^^ cm. The model is self-consistent 
for a wide range of Y as the saturation radius /?, = FTJq lies 
well within the photosphere with the radius Rp^ ^ 5 x 10'^ cm. 
The jet launching region should be located at several BH grav- 
itational radii. Assuming that the jet base radius is at Rq = 5Rg, 



,4v-3/2 



Mq, which 



we find the BH mass Mbh = (1.3-2.3) x lO^F, 
places the BH into the intermediate mass category. [Pe'er et aL[ 
(2007 1 argue that the ratio of energies is Yiq< 1. 



The reduced chi-squared x^/dof = 1.138 for the joint bkn- 
power fit and x^/ dof = 1 . 1 15 for the joint compPS fit are quite 
low despite good photon statistics. Even better fits down to 
X^/dof =1.0 can be achieved, when the host galaxy hydrogen 
co lumn density is allowed to vary between the time slices as 
in |Butler[ ( 2007| l. However, we do not see a physical reason 
for Nh to vary on the timescale of 10^ s, which corresponds to 
3 X 10'^ cm distance, while the characteristic star-formation 
size o f the host galaxy is R^q = 0.55kpc = 10^' cm (Svensso^ 
[et al.[ [2010). To test for potential discrepancies between the 
fits and the observations, we depict on Figure[T]the parts of fits 
and the normalized residuals for the 4-th time slice character- 
ized by the highest soft X-ray flux. The left panel shows the 
results for the bknpower model, while the right panel shows 
the results for the blackbody spectrum Comptonized by ther- 
mal electrons. The residuals do not systematically deviate 
from zero except slightly at the very low energies or for cer- 
tain absorption features. At the low energies Swift XRT re- 
sponse might be uncertain or the line physics might alter the 
spectrum. The most prominent absorption feature is around 
7.3 keV. 

4.3. Prompt Optical/UV Emission 

The prompt optical/UV emission observed by the UVOT 
instrument onboard Swift satellite was a matter of major dis- 
agreements between the former theoretical models. The same 
blackbody comp onent cannot produce bo th the optical and the 
X-ray emission ( Ghisellini et al.|2007b[ ), which is consistent 
with the physical picture within the tidal disruption scenario. 
For an estimate we take the magnetic field in the photosphere 
to be in equipartition with the jet kinetic energy. Then the 
Poynting energy flux is about Y /2 ^ 5 times larger than the 
total X-ray radiation flux. Hence the magnetic field in the 
photosphere has the strength 



B 



Ph' 



lO^G. 



(12) 



The compPS fit to the prompt X-ray spectrum shows that 
the non-thermal particles with the energies ^ 200 keV ex- 
ist in the photosphere. The fact that the emission contin- 
ues to the hard X-rays in BAT band reveals the existence 
of much more energetic particles with 7 ^ 10- 100, which 
readily emit synchrotron in the optical band. The synchrotron 
emission by an electron with a random Lorentz factor 7 30 
peaks at 300 nm for 10^ G magnetic field. Since the effective 
synchrotron cross-section is much higher than the Thomp- 
son scattering cross-section, the optical photosphere is much 
above the X-ray photosphere. In fact, very large jet parti- 
cle density makes the optical photosphere coincide with the 
front of the jet propagating with a bulk Lorentz factor F ^ 3 
into the surrounding medium. The jet front can propagate 
up to /?front = fobsr^c, wherc fobs is the time since BAT trig- 
ger. Continuous energy dissipation and particle acceleration 
at the front ensures there are enough energetic electrons to 
emit optical/UV synchrotron. Note that the interactions with 
the interstellar medium (ISM) at the early times, which lead 
to a forward shock, may contribute less to the emission. Low 
magnetic field at the jet front 

Bfront-IO^'G (13) 

at fobs = 3 X lO-' s is compensated by higher Lorentz factors 
7 ~ 300 of the electrons producing the optical/UV emission. 

The idea of the jet front emission explains the steep iv oc i'^ 
spectrum and the achromatic flux growth at different UVOT 
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TABLE 1 

Time-resolved spectroscopy of GRB060218 soft X-ray spectrum. Fitting with a broken power-law. 



Number 


Time period / [.s] 


Break energy 


High energy slope 


Absorbed flux 


Unabsorbed source flux 






£b,cak[keV] 


r2 


fabs[10"'erg s^'cm"-] 


^uiiabsllO^^erg s^'cm""] 


I 


164-478 


0.795 


1. 22b 


4.753 


5.234 


2 


478-691 


0.740 


1. 171 


7.411 


8.148 


3 


691-875 


0.820 


1.234 


8.242 


9.064 


4 


875-1049 


0.746 


1.290 


8.572 


9.575 


5 


1049-1226 


0.813 


1.419 


7.595 


8.575 


6 


1226-1414 


0.800 


1.519 


6.780 


7.795 


7 


1414-1620 


0.761 


1.650 


5.583 


6.631 


8 


1620-1854 


0.749 


1.816 


4.421 


5.467 


9 


1854-2119 


0.678 


1.962 


3.621 


4.783 


10 


21 19-2404 


0.712 


2.190 


2.940 


4.110 


11 


2404-2756 


0.656 


2.332 


2.265 


3.457 



Note. — The parameters of XSPEC hknpower model are shown. The low energy slope of the count spectrum is fixed at F i = - 1 . which corresponds to the Raleigh-Jeans tail of the 
blackbody emission with F,^ oc u' . The host galaxy hydrogen column density Nfj = 0.750 X lO^^cm"" with the met alljcity ^)7_otlsolar i Wiersema et al.|2007^ provides the best joint 
fit to all time slices with x~ = 5463.0 for dof = 4836. The galactic column is fixed at Nm = 9.4 x lO^'^cm"" followingjcampana et al.j^OO^ 



TABLE 2 

Time-resolved spectroscopy of GRB060218 soft X-ray spectrum. Fitting with the blackbody spectrum Comptonized by thermal 

ELECTRONS. 



Number 


Time period 


Photon temper- 


Electron 


Absorbed flux 


Unabsorbed 


Blackbody 


Lorentz factor 


jet base radius 




f [s] 


ature To[keV] 


temperature 


FabsLlO^'ergs^'cm"- 


source flux 


source flux 










7;[keV] 




Fiinabs [ 1 O^^erg s"' cm 


"^fsB[10-'ergs-'cm-- 




1 


164-478 


0.1039 


262.3 


4.586 


6.61 


2.630 


2.667 


0.906 


2 


478-691 


0.1052 


297.6 


7.159 


10.17 


4.033 


2.832 


1.091 


3 


691-875 


0.0994 


258.8 


7.965 


11.49 


4.465 


2.803 


1.245 


4 


875-1049 


0.1058 


230.2 


8.353 


12.29 


5.020 


2.898 


1.257 


5 


1049-1226 


0.0967 


176.6 


7.424 


11.52 


4.763 


2.744 


1.493 


6 


1226-1414 


0.1008 


145.8 


6.636 


10.84 


4.848 


2.753 


1 .560 


7 


1414-1620 


0.1119 


119.5 


5.517 


9.56 


4.784 


2.816 


1.487 


8 


1620-1854 


0.1116 


93.1 


4.369 


8.43 


4.686 


2.732 


1.733 


9 


1854-2119 


0.1324 


80.9 


3.614 


7.16 


4.412 


2.879 


1.393 


10 


2119-2404 


0.1489 


65.9 


2.970 


6.25 


4.284 


2.961 


1.275 


11 


2404-2756 


0.1421 


55.3 


2.268 


5.41 


3.963 


2.818 


1.485 



Note. — The parameters of the blackbody spectrum Comptonized by thermal electrons within a compPS model are shown. The optical depth is fixed at r^- = 1, which coiTesponds 
to Compton scattering at the photosphere. We find the host galaxy hydrogen column density Nh = 1.090 X 10^"cm~" in a self-consistent ioint fit over all time slices with x~ = 5358.9 
for dof = 4836 degrees of freedom. The Galactic column is fixed at = 9.4 X 10"°cm~" with solar metallicity following ^ampana et al. 1 2006 1. The quantity Yiq is the ratio of the 
to tal fireball energ y to the energy emitted in the X-rays in the units of 10. The last two columns are computed for Yiq = 1 or 7 = 10 within the thermal fireball emission model described 
in Pe'er et al.l|2007|. The jet base radius Rq ~ 10"^cm conesponds to the BH mass Mbh ~ 10"* Mq assuming the jet is launched within several from the center. 




Fig . 1 . — The parts of XSPEC spectral fits and the normalized residuals for the 4-th time slice. This time slice corresponds to the highest soft X-ray flux. The 
left panel shows the results for the bknpower model, the right panel shows the results for the blackbody spectrum CoiTiptonized by theiTnal electrons within the 
compPS model. The model parameters can be found in Tables0and[2] 
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frequencies till ^ 3 x 10'* s. The jet front is transparent to the 
hard UV and X-ray photons. Since the particles are subject 
to rapid cooling, the spectrum of re-accelerated particles may 
substantially deviate from a power-law and resemble a very 
hot Maxwellian. Then the resultant self-absorbed spectrum 
could be the Rayleigh-Jeans part of the blackbody oc 
as opposed to the absorbe d synchrotron oc v^l'^ spectrum 
( Rybicki & Lightman 1979 ). The obs erved te mperature of the 
optical/UV radiation is about 10'' K ( jGhiselUni et al. 2007a1), 
which violates the energetics of the source in a non-relativistic 
emission model with the total emitted energy Ebb ^ 10''' erg. 
In turn, the relativistic jet front emission does not violate the 
source energetics. Relativistic Doppler boosting with F ^ 3 
reduces the emission temperature down to < 3 x 10^ K, which 
readily leads to Ebb ^ 10^" erg. The inferred large optical 
emission radius 10''*"'^ cm at 3 x lO-' s is consistent with Rf^vA 
for the mildly relativistic expanding shell. A more detailed jet 
model for optical/UV emission by Ghisellini et al.| ( j2007b| ) is 
similarly consistent with the source energetics. 

4.4. Temporal Analysis 

The spectral modeling allows for proper determination of 
fluxes at different times to study the temporal behavior of 
the source. Figure |2] shows the absorbed and the unabsorbed 
source fluxes as functions of time. The soft X-ray flux peaks 
at around 800 s following Swift BAT trigger We fit the 
lightcurve with the scaled accretion rate dependence on time 
M(f) for a tidal disruption of a MS sta r by a SMBH Such a 
dependence for /Jj- = 1 was taken from |Laguna et al!] ( 1993| l. 
The best fit corresponds to a time delay of Af = 1810 s be- 
tween the tidal disruption and Swift BAT trigger. The peak at 
fpeak ~ 2600 s corresponds, for example, to the tidal disrup- 
tion of a m^D = OJSMq WD by a 1 x IO^Mq IMBH (or a 
mwD = OMMq WD by a 2 X IO^Mq IMBH) at I3t = 1. The 
tidal radius in this case is l6Rg (or 8/?g). A deeper encounter 
with (3t > I may be needed to initiate a nuclear burning for 
an accompanying supernova, but the dependence of critical 
f3x, for which the effective nuclear burning starts, on the WD 
and the BH masses is not known at present. The encounter 
in question may indeed have very small pericenter distance 
with /jj- ^ 1 and not violate the peak time c onstraint, as |De| 
IColle et al.| ( [20T2l ); [Guiilochon & Ramirez-Ru iz (2012| fiiid 
very weak dependence of peak times on /3t- 

The accretion rate curve M{t) provides a good fit to the 
unabsorbed source flux. However, such fitting procedure is 
prone to multiple caveats. First, the accretion disk onto the 
BH takes a finite amount of time fmagn ^ lOOfoit ^ 250 s to 
generate the strong poloidal magnetic field. Despite the mag- 
netic field generation time is much less than the characteristic 
dynamical time fmagn ^ fdyn ^ 1000 s, the precise fit is not ex- 
pected. Secondly, the transition from the accretion power Mc^ 
to the emission power in a certain ban d is not trivial. Th e jets 
are expected to be radiatively efficient ( |Pe'er et al.|2007| l, and 
most of the emission, especially at the late times, falls into the 
0.05-10 keV band. Only if the jet kinetic power is a constant 
fraction of the accretion power, then the soft X-ray luminos- 
ity is expected to approximately follow the accretion power. 
Finally, the GR effects are not included in these estimates. 
|Haas et ar] ( j2012j emphasized that the ultra-close encounters 
with small ~ several x /?„ may lead to non-trivial temporal 

behavior of M. 

The estimated large BH mass and the small Lorentz factor 
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Fig. 2. — Soft X-ray lightcurve of GRB060218 fitted witii the accretion 
rate temporal evolution M(t). Shown are the unabsorbed source lightcurve 
for the compPS model (dark solid) and for the bknpower model (red dashed), 
the observed absorbed lightcurve for the compPS model (lower dotted) and 
for the bknpo wer model (up per dotted), and the scaled fallback accretion rate 
based on Laguna et al.H199 3 1 (light solid). The scaled accretion rate is offset 
by At = IblU s, which the most bound debris takes to reach the BH. Two 
thin lines (t + Aty^ and {t + Aty^/^ are shown to emphasize the asymptotic 
behavior of the absorbed flux and the unabsorbed flux, respectively. While 
the observed flux decays as F oc r^, the unabsorbed flux exhibits F oc r^/^ 
late-time behavior in agreement with the tidal disruption scenario. 

may lead to a smooth lightcurve. The lightcurves of typical 
GRBs are highly variable, which is attributed to fast variations 
near small BHs with masses Mbh ^ IOM0, when the BHs 
launch shocks at different speeds. The variation timescale is 
correspondingly longer for the BHs with much larger masses 
Mbh ^ IQ'^Mq. Those variations amplify, when the shocks 
moving at different high Lorentz factors collide. The energy 
release and the radiated energy are relatively small, when the 
shocks with small F's collide. Lastly, the collisions of the 
slow-moving shocks are likely to happen inside the photo- 
sphere, so that any variations are smoothed. Thus, a slow 
moving jet of GRB060218 is expected to vary on a very long 
timescale and the variations are expected to be small. The 
object should exhibit a smooth lightcurve. On Figure [3] we 
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Fig. 3. — Log-smoothed periodograms I{T) jPapadakis & Lawrence] 1993^ 
of the lightcurves binned to 1 s. The flat red curves depict simulated white 
noise, GRB060218 is shown by a dark/blue curve, while 16 other GRBs are 
shown by light/green curves. The white noise has the same mean count rate 
and total counts as GRB060218. The other GRBs are selected to have the 
peak count rate between 142cts s"' and 166cts s"', while the peak count rate 
is 148cts s"' for GRB060218. The break time tt, at which the periodogram 
starts deviating substantially from the noise, corresponds to a characteristic 
system timescale. The break time for our event is ti, ~ 300 s, which is con- 
sistent with the emission rise timescale in the tidal disruption. Other GRBs 
have substantially lower break times = 5 — 40 s, and the distribution of ti, 
does not continue to 300 s. This indicates a different origin of GRB060218. 
Log-smoothing is performed to 0. 1 1 dex. 

present the log-smoothed periodograms 1(7) ( [Papadakis &| 
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Lawrence||T993] l of GRB060218 and other lightcurves. The 
dark (blue) rising curve is the periodogram of GRB060218 
prompt X-rays. The dark (red) constant curves are the peri- 
odograms of 20 white noise implementations with the same 
mean count rates and total counts as GRB060218. The light 
(green) curves represent the selection from the online Swift 
catalogue of 16 GRBs, whose peak photon count rates are 
the closest to the peak count rate 148cts s"' of GRB060218. 
The white noise implementations are flat as expected, while 
our candidate tidal disruption source rises above the noise at 
a very long break time of about f/, = 300 s. Other 16 GRBs 
show much smaller break times of 5-40 s. For some of 
the sources the break time coincides with the characteristic 
flux decay time, while most sources exhibit variability on a 
shorter timescale. As predicted, the lightcurve of GRB060218 
is smooth and shows no variability faster than the emission 
rise timescale. This strengthens the case for the source being 
a slow dense jet from an IMBH. 

4.5. Steep Decay Phase and Afterglow 

The Swift satellite was unable to observe the source from 
2780 s till 5900 s counting from the BAT trigger. The ob- 
served X-ray flux at 5900 s is very low, which is seemingly 
inconsistent with t~^l^ behavior of the accretion rate. How- 
ever, the blackbody source flux may still be large, while the 
observed flux may be a factor of 40 lower due to the action 
of absorption and incomplete overlap with the XRT band. 
The XRT instrument observes mostly Comptonized photons 
with its effective energy range 0.3-10 keV. The energy of 
the Comptonizing electrons goes down with time during the 
prompt phase, which translates into the ratio of unabsorbed 
to absorbed fluxes rising from 1.4 to 2.4 by 11-th time slice 
(see Table [2]). The ratio of fluxes keeps dropping till lO'* s 
in the steep decline phas e. After that muc h shallower de- 
cline F oc f"'-^ follows (Camp ana et al.|[2006 ) with th e spec- 
trum consistent with a power-law F^, (x v~^'^ (Soderberg et al. 
|2006| l. While the temporal decay of this late phase is con- 
sistent with the afterglow, the spectrum is substantially softer 
than th e typical afterglow spectrum emitted by the external 
shock ( Toma et al.||2007[ l. The softer spectrum and the long 
afterglow may be explained by the late activity of the source 

toderberg et al. 2006 ). This is an especially viable idea, since 
e jet kinetic luminosity is expected to stay high for a long 
period of time following a tidal disruption. 

We model the steep decay spectrum within the time inter- 
val f = (5950,7070) s with a combination of the blackbody 
component Comptonized by thermal electrons and the power- 
law emission from the external shock. We fix the absorption 
at the level determined with compPS fitting of the prompt X- 
ray emission. The fit allows for a large range of power-law 
slopes without a significant change in x^. We fix the power- 
law slope at Tpi = 2 (/v oc t /~^), which corres ponds to the af- 
terglow of a typical GRB ( [Soderberg et al.||2006) . The ab- 
sorbed combination of the compPS and the powerla w mod- 



Butler 



els leads to = 43.1 for dof = 36. As suggested by ^ 

(j2007 ), the spectrum contains lines, which become prominent 
at the late times. Lines is the natural outcome of the atomic 
processes in a warm absorbing/emitting jet baryonic material. 
The jet baryonic material consists of WD debris rich in oxy- 
gen and carbon and practically devoid of hydrogen. The ab- 
sorption by cooled down jet exhaust may become strong at 
the late times. We model such absorption by a blueshifted 
zvphabs model, where we set the abundance of a single chem- 



ical element to 1000 and the abundances of the other elements 
to zero. The model with a pure oxygen cold absorber con- 
verges to a blueshift z = -0.422, which corresponds to the 
bulk Lorentz factor T = 1.73 consistently with the fits to the 
prompt spectrum. The fit with = 37.3 for dof = 36 and 
the normalized residuals are depicted on Figure [4] The fit 
can be marginally improved with the addition or iron ele- 
ments. We do not present such modeling here for simplic- 
ity. The addition of up to 30% of carbon does not change the 
fit. The best-fitting oxygen column density A^o = 1 x 10"* cm"^ 



dalit and I'nEdcd model 



^ 0.2 - 




Fig. 4. — Best fit to the spectrum at the late times t = (5950,7070) s with 
= 37.3 for dof = 36. The model consists of the blackbody emission Comp- 
tonized by thermal electrons absorbed by blueshifted oxygen and a powerlaw 
with a fixed slope F,,/ = 2 (iv oc v^^). The sum of these two components is 
absorbed by the galactic Nh = 9.4 X lO^^cm"^ and Nh = 1.090 X 10^-cm"^ 
of the host galaxy with their respectively metallicities all fixed from the 
the early time modeling. The blueshift z = —0.422 corresponds to the bulk 
Lorentz factor F = 1.73. The best-fitting cold oxygen column density is 
No = 1 X 10'**cm"^, which can readily be provided by a cooling jet exhaust. 
The blueshifted oxygen manifests as an absorption feature around 1 keV.The 
unabsorbed source flux in f|inabs(0 05 — lOkeV) = 1.8 X Kr'erg s^'cm"^, 
which coincides with the extrapolation based on t^^l^ law, while the heav- 
ily absorbed observed flux is Fj^,^ = 4.0 X 10""erg s^'cm"^. 



can be readily provided by the cooling jet material. Very soft 
blackbody has a temperature T = 0.052 keV. The best-fitting 
flux is Funabs(0.05- lOkeV) = 1 .8 x lO-'^erg s-'cm'^, which di- 
rectly traces t~^l^ decay law from the time of peak flux, while 
the observed absorbed flux Fabs = 4.0 x 10~"erg s"'cm"- is 
40 times lower. The absorbed power-law contribution to the 
total source flux is Fabs.pi = 2.5 x 10~''erg s"'cm"^, which 
is consistent with the emission at later times being domi- 
nated by the po wer-law. T he application of the photospheric 
emission model ( |Pe'er et al. 2007) to determine the Lorentz 

factor and the jet base radius gives r = 1.45710'''* and Rq = 
1.1 X lO'^F"-'/^ cm, respectively. A larger ratio of total to 
radiated energy F = 4 x Fprompt makes the numbers consistent 
with the results from the prompt phase. Since the blackbody 
is too soft to be directly observed, there is a substantial de- 
generacy between its temperature and normalization, which 
leads to the unabsorbed flux being uncertain by a factor of 3. 
Taking this uncertainty into the account, the unabsorbed flux 
^iinabs is Consistent with the t~^l^ law even for larger Y . We 
might not be able to compute a more self-consistent physi- 
cal model of the steep decay phase unless we incorporate the 
atomic physics of warm emitters/absorbers, which lies beyond 
the scope of the present paper. 
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According to Cannizzo et al. (2011 1 switching from t 
fallback dominated behavior to C^^^ disk spreading behavior 
of the accretion rate occurs at 



f,= 1.25xl0"'s 



-3/2 



Rq 



3/2 



Mn 



(lO^Mg 



{ — ){- 
Vo.l/ V2/3' 

(14) 

where = Rdrc/Rp is the effective extent of the disk. Tak- 
ing a WD with a mass O.86(O.75)M0 and a radius R^, = 
7.1(8.1) X 10*^ cm disrupted by a 2.0(1.0) x lO'^M© BH we 

estimate the transition time to be = 2.6(1.6) x 10"* /Sj-^''^ s 
for 77 = 2. Thus, the transition to the shallow disk spreading 
behavior occurs around the transition to the afterglow phase. 
The temporal index of the power-law flux decay -1.2 ±0.1 
during the afterglow ( Soderberg et al. 2006 ) is consistent with 
-1.33 index of the mass accretion rate at the late times. 

4.6. Associated Supernova and Host Galaxy 

Another challenge to the WD/IMBH tidal disruption model 
is the associated supernova. As discussed in § |2| the super- 
nova ejecta mass Mq ^ (1 -2)M0 is consistent with the high 
end of the WD mass distribution. The unbound mass fraction 
is influenced by the interplay between the energy spread dur- 
ing the tidal disruption and the energy release in the supernova 
explosion. The outflow velocity of tidal disruption debris is 
20 X lO^km s"' for the WD and the BH masses estimated 
above, which is comparable to the typical outflow velocity in 



a Type la supernova (10-20) x 10 km s (Wang et al 



2009 I. 



Thus, the supernova SN2006aj with the en ergy release corn- 
parable to that of a Type la supernova ( Kho khlov et al.|1993) 
can unbind most of the WD material, which then contributes 
to the supernova ejecta mass. Only a small fraction of debris 
accretes back onto the BH. The nebular phase of SN2006aj 
exhibits strong oxygen emission (Mazzali et al. 2007), which 
is consistent with the ignition of a carbon/oxygen WD. 

A supernova explosion accompanying a tidal disruption is 
expected to be heavily asymmetric. In particular, there could 
be a velocity shift of supernova lines with respect to the ve- 
locities of the host galaxy. Yet, the velocities and the re dshift 
z ~ 0.033 were only estimated for the host galaxy (Mirabal' 
et al.|2006{ Modjaz et al.|2006| ) and no measurements of su- 
pernova velocities exist. 

Since the IMBHs are expected to reside in the nuclei of 
dwarf galaxies, the other major test is the coincidence of the 
supernova position with the center of the host galaxy. A 
star-formation radius for GRB060218 host is /J^o = 0.55 kpc 
(Svensson et al. 2010), which corresponds to a half-light ra- 
dius R50 « 0.45 kpc. The supernova and the host galaxy were 
observed with Hubble ACS instrument at different epochs. 
Figure|5]shows the position of the supernova (green contours) 
relative to the intensity plot of the host galaxy at a late time. 
The host galaxy has a somewhat ir regular morphology as can 
be seen on images by Misra et al. (201 1 ), which complicates 
finding its center. The uncertainty of center determination is 
about dR < O.larcsec « 70 pc. The supernova position ap- 
pears to be consistent with the center to within O.larcsec. 
The ratio of the uncertainty dR to the half-light radius is 
dR/R^o < 0.15, which is much less than a normalized off- 
set dR/Rso = 0.98 fo und for a representative set of long GRBs 



( [Bloom et al.|2002 



It is very rare for a supernova to happen 

that close to a galactic center at random. 

The stellar m ass of the host gal axy is estimated to be 
Mst ~ 10^ ^M© (Ferrero et al. 2007 l Observations suggest 
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Fig. 5. — Position of SN2006aj supernova (green contours) witliin the 
host galaxy (grayscale intensity plot) from Hubble ACS observations in red 
F814W band. The supernova position is consistent with the center of the 
host. The supernova image is computed as the subtraction of the late time 
(26 Nov 2006) image from the early time (18 Jul 2006) image. The late time 
image without supernova subtraction is taken to represent the host galaxy, 
which has in'egular morphology. The subtraction is not performed to avoid 
the oversubtraction of the galactic nucleus were the supernova to happen at 
the very center. 

that the bulge mass to the total stel lar mass ratio (B/T) of 
galax ies is within 15- 100% range ([Schramm & SilvermanI 

'Wa 



2012| l and B/T = 15% for the Milky "Way (McMillan 201 IT 
I'hen, extrapolat ing Mbuige ~ Mbh relation ( jMagorrian et ak] 
[19981 [Marconi & Hunt|2003| l 



Mbh/M^,,^, ~ 0.002, 



(15) 



to low masses, we get the central BH mass to be Mbh ^ 
(0.5-3) X 1O'*M0. This estimate, though much less certain, is 
consistent with two previous estimates based on the jet ther- 
mal emission and the lightcurve peak time. 

4.7. Accretion Flow and Jet Energetics 

Let us compute the energetics of the jet for the estimated 
BH mass. Assuming the supernova does not unbind the gas, 
the accretion rate onto the BH peaks at 

M « 2 X lO^Meyr"' , (16) 

which corresponds to the accretion power Mc^ « 1 x 
10^"erg s"'. Then the Blan dford-Z najek jet power is Pbz ^ 
O.OlMc^ - 1 X lO'^'^erg s"' (j McKinney 2005 ) for the low ef- 
ficiency associated with simulations initiated with the weak 
magnetic field. The isotropic equivalent radiation power for 
the opening angle 9 = 0.5 and the ratio F = 10 is ^'rad.iso ~ 
2 X 10^'^erg s"'. This is substantially higher than the observed 
peak isotropic luminosity Pobs.iso = 3 x 10'*^erg s"'. The dis- 
crepancy can be partially explained by a small fraction of a 
WD material falling back onto the BH, since the supernova 
may unbind most of the gas. Other potential explanations in- 
clude GRB060218 happening substantially off-axis, the WD 
being mostly unb ound in an encounter with a spinning BH 
^Haas et al.|2012| l, or the source failing to launch a powerful 
jet. The latter is especially viable, since the source has a very 
limited amount of time to generate the ordered magnetic field. 
A poloidal magnetic field 10 times weaker than the equipar- 
tition value given by the equation ([8]l leads to consistent jet 
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energetics. The disk may launch a mildly relati vistic outflow 
by itself via the Blandford-Payne mechanism (Blandford & 



Payne 1982) i. An off-axis GRB is inconsistent with the ab- 
sence of a jet break (fSoderberg et al. 2006 ). 

The X-ray emission is dominated by the super-Eddington 
jet at the early times, but the accretion disk may be brighter at 
the late times. Our disruption models have the peak accretion 
rate of about Mpeak ~ 2 x lO^'Moyr"' at the time 2600 s after 
the disruption (neglecting the effect of blowing the material 
away by a supernova). Then assuming t~^l'^ law we find that 
the accretion rate is down to the Eddington value for a stan- 
dard e = 0.1 efficiency at fedd ^ 1 yr. As described in § 4.5 
the accretion rate is expected to follow a shallow f""*/^ slope 
after about fj; ^ 2 x lO'* s due to the viscous evolution in the 
disk. Yet, as we will describe below in § |5.1[ when the outer 
disk boundary diffuses out to be much larger than ^ 20/?g, 
then the energy transfer between the inner and the outer disk 
kicks in. The energy from the inner flow unbinds the outer 
material and the acc retion rate onto the BH dramatically de- 
creases. Cannizzo & Gehrels (2009) derived the temporal de- 
pendence of the accretion rate in the RIAF with the energy 
transfer/convection to be f"4(i+*)/3^ where the parameter s con- 
trols the convection efficiency. The value of s varies from 
i = for no energy transfer up to i = 1 for the maximum en- 
ergy transfer power Then the accretion rate may decrease as 
steeply as t~^l^ in the RIAF phase. We w ill discuss the RlAFs 
with convection in more detail in § |5.1| The onset of RIAF 
with the energy transfer may explain the steeper XRT flux de- 
pendence F (X f"'-^^ at times f > 2 x 10^ s, accompanied by 
a nearly constant hardness ratio. Therefore, the accretion rate 
may fall below the Eddington rate at times f ^ 1 yr 

Yet, the Eddington luminosity for a 2 x IQ'^Mq BH cor- 
responds to a flux Fx ^ 10~'^erg s"'cm"^, which is about a 
factor of 100 higher than the late-time flux observed by the 
XRT. This is not a contradiction, however, as at times much 
before 1 yr the accretion disk might not be visible at all due to 
obscuration by the outflowing debris. As discussed in our pre- 
vious paper ( [Haas et al.|2012| l, the outflowing debris consist- 
ing of carbon and oxygen have high absorption cross-section 
cr « 2 X lO'^^cm"^ across 0.5- 10 keV band. The obs erved 



matter velocity in SN2006aj is about 1.5 x lO'^km s (Maz- 



zah et al. 2006b). A tidal disruption may scatter matter in par- 
ticular directions leaving huge voids, though which the inner 
disk could be observed. A tidal disruption of an object with 
the size comparable to the tidal radius rather tends to sc atter 
the debris over the large solid angle (| Haas et ar]|2012[ ). In 
addition, the supernova, whose energy release is comparable 
to the energy release in a tidal disruption, should sca tter the 
debris more uniformly in all directions. Following Haas et al. 
( |2012 i we find that the absorption optical depth unity r = 1 
is achieved around t ^ 3 yrs, at which time the source is ex- 
pected to rebrighten. As the accretion disk should have a sub- 
stantially sub-Eddington accretion rate by 3 yrs, the emitted 
flux levels Fx < 10~'^erg s"'cm"^. The accretion flow with the 



sub-Eddington a ccretion rate settles into a thin disk (Shakura 
& Sunyaev|1973) l, which has the inner temperature about 



GMbhM 



1/4 



(17) 



The temperature is Tin ~ 0.15 keV for the accretion rate 
equal to 10% of the Eddington value and the BH mass 



2 X 1O'*M0. Very soft blackbody spectrum absorbed by the 
galactic column and the host galaxy column leads to the ratio 
^unabs/^abs ~ 10 of the modcl flux to the observed absorbed 
flux and the observed flux F^^s ^ 10~'"^erg s"'cm"^, which 
might not be detectable. A tantalizing possibility that the disk 
dominates the emission starting as early as 10** s, at which 
point the source flux equals the Eddington flux of a hypoth- 
esized IMBH, is unrealistic. First, the optical depth of su- 
pernova debris is about r ^ 10'* at that time, and the special 
directions of low t < 1 are not expected. Second, the accre- 
tion rate at 10"^ s is expected to be highly super-Eddington, if 
only a tiny fraction of the WD material is bound to the BH. 
Third, even if the accretion rate at 10"* s was at the Eddington 
level, then the extrapolated jet power at 10^ s for the standard 
jet efficiency would be ^ 100 times lower than the true soft 
X-ray power for a jet with 9 = 0.5 opening angle. 

5. EVENT RATES 

The IMBHs are thought to reside either in the GCs or in 
the dwarf ga laxies. As we discussed in the earlier work (Haas] 
|et al.||2012"] l, the d isruption rate of stars by IMBHs in GCs 
IS very uncertain. [Baumgardt et al.| ( |2004| l predicts that for 
a GC with a IO^Mq central IMBH, the optimistic disruption 
rate of stars is 10~^yr"' per GC, while 15% of all disruptions 
are those of WDs. Then, following McLaughlin (jT999]l we 
theoretically estimate the space density of the GCs to be 

nccth « 34Mpc-^ (18) 

which leads to the rate ( |Haas et al.||2012| l Rwd-imbh ^ 
500yr-'Gpc"^ of the WD disruptions by flie IMBHs. Obser- 
vational constraints on the GC populati on suggest a slightly 
low er space densit y of the GC s (Brodie & Strader] |2006[ 
|Ramirez-Ruiz & Rosswog|2009| l 



nGC,obs~4Mpc , 
so that the rate estimate should be revised down to 



Rn 



■50yr-'Gpc"l 



(19) 



(20) 



The disruption rate of stars in dwarf gala xies was esti- 
mated to be very high. If M-g (IGebhardt e t al. 200()|l_an_d 
Mbh -A^buige (|Magon-ian et al.|1998', 'Marconi & Hunt 2003 



1 



relations hold down to the low BH masses, then the predictei 
tidal disruption rate is hig her for the dwar f galaxies compared 
to the other galaxies (Wang & Merritt|2004) . Dependent on 
the radial stellar profile the disruption rate in a dwarf galaxy 
with a 1O'*M0 BH varies from Rimbh = 10~^yr"' per galaxy up 
to Rimbh = 10~'yi~' pei" galaxy. For an estimate we take the 
rate of Rimbh = 10~''yr"' per galaxy. Note that during a tidal 
disruption some fraction / of a WD accretes onto a BH. Then 
the mass of an IMBH grows at a rate 



dMsH 
dt 



■ fRlMBHniwD- 



(21) 



Taking a large fiducial fraction / = 0.5 we find that the 
BH mass doubles every tj = 2 x 10^ years solely due 
to tidal disruptions of WDs, where we assumed that the 
same Rwd-imbh /Rimbh = 15% of the disrupted stars are the 
WDs. This doubling is consistent with the estimated a ge of 
GRB060218 host galaxy 200 Myr (Ferrero et aL][2007j i. As 
we will show in the next subsection, an IMBH swallows only 
a small fraction of a MS star, so that the IMBH mass growth 
rate due to the MS stars may not be dominant. 
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According to observations ( {Ferguson & Sandage|1991 1 the 
number density of dwarf d galaxies with V magnitude My > 
-15.5 is about the number density of large elliptical E and 
spiral S galaxies. Thus, the space density of dwarf galaxies is 
"dwarf ^ 10~^Mpc~^. Then, assuming that every dwarf galaxy 
has a small enough IMBH, we arrive at a rate estimate 



1500yr-'Gpc-^ 



(22) 



which is much larger than the rate in GCs. One event per year 
for the field of view (FOV) 10% of the sky would happen 
as close as lOOMpc. Such disruptions should be detected by 
Swift BAT instrument with FOV of approximately 10% of the 
sky. Collimation into a wide jetted outflow may reduce the 
rate down to one event per several years consistently with a 
single GRB060218 source at 143 Mpc distance. 

Since the WD disruptions by the IMBHs produce super- 
novae, we c an comp are the rate of those with the total super- 
nova rates. |Li et al.] ( |20Tl) find the local rate of supernovae 
Type la to be 



/?s^v,/a«3x 10V"'Gpc-^ 



(23) 



and the rate of supernovae Type Ib/c to be RsN.ibc ~ 3 x 
10"^yr"'Gpc~"'. In an optimistic estimate, if a supernova is pro- 
duced in every tidal disruption of a WD, then one disruption- 
induced supernova happens per 20 Type la supernovae. A 
distinctive feature of disruption-induced supernovae is their 
location in the galactic nuclei substantially close to the photo- 
metric centers of their dwarf host galaxies. This rate equals to 
the rate of calcium-rich gap transients, whose properties may 
be explained by the u nderluminous explosions of the WDs 
(Kasliwal et al. |2012| l. However, those transients occur very 
far from the galactic nuclei. 

5.1. Disruptions of MS Stars by IMBHs 

In the present paper we consider the radiative signatures of 
the WD disruptions by the IMBHs. However, the disruptions 
rates of the MS stars by the IMBHs are several times higher. 
Then the question is whether these more frequent events pro- 
duce equally distinct and observable signatures. In this sub- 
section we show that MS star/lMBH disruptions are quite dif- 
ferent: they are much longer and much fainter Despite lower 
rates, it might be easier to observe the WD/IMBH disruptions 
compared to the MS star/lMBH disruptions. 

The fallback timescale is about ffb ^ 20 days for a disrup- 
tion of a IMq MS star with a solar radius by a IO'^Mq IMBH 
with a pericenter radius equal to the tidal radius Rp=Rj-- The 
peak fallback rate is Mpeak ^ 5M0yr~'. While a tidal disrup- 
tion radius is Rt ^ 2QR„ for a WD, a MS star with the solar 
mass and the solar radius has a much larger Rj ^ 1000/?g. 
Since the peak fallback accretion rate Mpeak ^ 5M0yr"' is 
much above the Eddington rate MEdd ^ 2 x lQ~*MQyi~\ the 
RIAF settles from the outer disk radius 7?disk ~ ^Rp down to 
the BH. A RIAF with an outer radius greater than about 20Rg 
settles into the convection dominated accretion flo w (CDAF) 
( [Narayan et al.|2000[[Quataert & Gruzinov|20"00| , where the 
energy transport between the inner and the outer flow starts 
to play a role. Unlike the fallback disks following the WD 
disruptions, the disks following the disruptions of the MS 
stars always exist in CDAF state. A flow with convection has 
a shallow d ensity profile p oc down to the inner radius 
i^in ~ 20Rg ( Abramowicz et al.|2002^ . Only a small fraction 



of available matter accretes in CDAF state. The sustained BH 
accretion rate is 



M = Mfb(-RdiskMn) 



1.5-13 



(24) 



The density slope /3 = 1 .5 -s = 0.5 - 1 .0 was found in the 
numerical simulations (see [Yuan et al. 
The va lue f3 = 0.8-0.9 was estimated 
( |2012p for Sgr A*. The flow settles on a viscous timescale 
fvisn ^ 1 day, which is shorter than ffb ~' 20 days. Thus, the 



2012' for t he review). 
ly Shcherbakov et al. 



peak accretion rate of the CDAF is given by the formula ( 24 1 
with Mfb = Mfb.peak- The peak accretion rate onto the BH can 
be a factor of 10- 100 lower than the peak fallback rate 



Mpeak- (0.01 -0.1)Mft .peak- 



(25) 



Let us make a fiducial estimate for a BH mass 10'^Mq and a 
disruption of a Sun-like star for Rp = Rj assuming density 
slope (3 = 0.85. In this case the accretion rate is lowered 
by a factor of « 23 due to the action of convection accord- 



ing to the formula (24i, and the resultant jet power is Lkin = 
O.OlMpeakC^ - lO^'^erg s"' with the same fiducial efficiency 
Pq ^ 0.01 used in estimates for the WD disruptions. A jet with 
such low power has the photospheric radius smaller than the 
saturation radius. The radiation decouples from matter before 
the jet can accelerate, and the powerful photospheric emission 
is not expected. The photospheric radius corresponds to the 
bulk Lorentz factor F « 2. If the source is beamed towards 
us, then the observed temperature is Tob ~ 0.5 keV. The ra- 
diated blackbody luminosity is Lbb — 1 x 10'*-' erg s"', which 
corresponds to the distance d ^ \Q Mpc to match the peak 
flux level Fpeak ~ 10"** erg s'^cm'^ of GRB060218. The ab- 
sorption by the cooled down jet exhaust or the dense wind, 
which consist mostly of hydrogen and helium, may substan- 
tially lower the X-ray luminosity and turn the source into an 
optical/ UV trans ient similar to the disruption of a MS star by a 
SMBH dS trubbe & Quataert 2009; Str ubbe & Quataert|201 \) . 
The optical/UV transients lasting for 20 days and releasing up 
to lO'^^erg s"' may be confused w ith variability of the AGNs. 
According to the stand ard theory ( Cannizzo & Gehrels|2009[ 
Cannizzo et al. 112011) at t^ ^ 100 days the source should 



switch from the fallback-dominated f"^/^ behavior to the disk 
spreading behavior f^C+^V^ f-2.2 yet, the timescale fv was 
derived for the adiabatic RIAFs. The temporal behavior of 
radiatively inefficient fallback accretion disks with convec- 
tion were not explored to our knowledge and should be in- 
vestigated in the future work. If absorption is inefficient, then 
one transient per year is expected at a distance li 100 Mpc 
with the peak flux level of about Fpeak ~ 10~"*erg s"'cm~^. 
Searches within the existing and future X-ray and optical 
surveys should identify such sources and constrain the rates 
of tidal disruptions of various kinds. ROSAT All-sky sur- 
vey search revealed several nearby tidal disruption candidates 
with redshifts z < 0.1 (Komossa,2002^ . Some of those candi- 
dates, e.g. RX J1242-ni9, cannot be easily matched with a 



large host galaxy ('Komossa & Greiner|1999 Komossa et al. 
2004| l. The absence of suitable X-ray/optical transients in the 
data might not immediately invalidate high tidal disruption 
rates. For example, if the jet power in the disruptions of the 
MS stars is as weak as in GRB060218, then the X-ray/optical 
flux from such disruptions is expected to be up to 100 times 
lower. The detailed analysis of the tidal disruptions of the MS 
stars is beyond the scope of the present manuscript. 
A candidate source for a stellar disruption by an IMBH is a 
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source in a GC NGC1399 explored by |Irwin et aT]d20T0| . Ob- 
served at a distance 20 Mpc the source exhibited X-ray lumi- 
nosity of 2 x 10-''' erg s"' and luminosity of a few x 10^^ erg s"' 
in oxygen and nitrogen optical lines. These properties were 
found consistent with a disruption of a hor izontal branch star 
by a 50- lOOM© BH ( jClausen et al.|2012| i. 



6. DISCUSSION AND CONCLUSIONS 

In this paper we explore the theory of the WD tidal dis- 
ruptions by the IMBHs, an understudied, but a very promis- 
ing type of an encounter. The high rate of the WD tidal dis- 
ruptions in dwarf galaxies warrants the search for candidates 
among already observed objects. We identify GRB060218 
with the accompanying supernova SN2006aj as a promising 
candidate and model its temporal/spectral properties and the 
properties of the host galaxy. GRB060218 was explored pre- 
viously within the supernova shock breakout model and the 
model, where a jet is launched by a newborn BH or a neutron 
star. We find the tidal disruption model to be a viable alter- 
native, which more naturally explains some features of the 
candidate. In this section we briefly discuss the application of 
the tidal disruption model to GRB060218/SN2006aj source, 
and compare different models of the source. 

The tidal disruption model performs equally well compared 
to other models in terms of event rates. The appearance of 
GRB060218 at a redshift z = 0.033 is consistent with a dis- 
ruption rate 10~^yr"' of all stars in a dwarf galaxy. The rate of 
low-luminosity GRBs such as our sourc e was estimated to be 
about 10% of supernova Type Ib/c rate ( Toma et al.|2007 1. A 
number of shock breakouts were observed to accornp any su- 
pernovae at all distances ( Woosley et al. 1999; Mazza li et al. 
|2006a| [Chevalier &"F ransson 2008 ; Schawinski et al. 2008P , 
and the sample of shock breakouts is likely incomplete. 

The long emission duration 2600 s provides a mass estimate 
1 X IO'^Mq for an IMBH within a tidal disruption scenario. 
The relat ivistic shock breakout model also explains the event 
duration (Nakar & Sari|2012| i. The extreme smoothness of the 
spectrum is consistent with either the supernova shock break- 
out or the jet launched by an IMBH. The jet launched by a 
stellar mass central object may produce more time variability. 

While the observed absorbed X-ray lightcurve of the event 
can be fitte d with the exponential d ecay within a shock break- 
out model (ICampana et al.p006[ l, the source flux corrected 
for the absorption exhibits (f - fo)~^^^ behavior specific to the 
tidal disruptions. Moreover, the full soft X-ray lightcurve in- 
cluding the rise and the decay phases can be fitted well with 
the accretion rate temporal dependence M{t) for a tidal dis- 
ruption. The steep flux decay at ^ 6500 s can be ascribed to 
the source being more heavily absorbed at late times, while 
the jet power and the radiation power still scale as f~^/^. The 
power-law decay is not expected in a shock breakout, while a 
shallower tem poral slope is expected for a jet-powered source 
( [Soderberget al. 2006 ). The temporal slope of the afterglow 
phase, powered by the central engine, is consistent with the 
accretion rate temporal dependence r"*^^ in a RIAF dominated 
by the disk evolution behavior. 

The observed blackbody flux and the blackbody temper- 



ature can be readily modeled by the jet photospheric emis- 
sion to provide the estimates for the IMBH mass 1 x IQ'^Mq 
and the Lorentz factor F ^ 3. Moreover, these estimates stay 
consistent between different time slices in our time-resolved 
spectroscopic analysis. The rise of the blackbody temperature 
accompanied by the constant blackbody flux and the dropping 
total flux leads to the same BH mass and Lorentz factor over 
many time slices within the fireball model. The relativistic 
shock in the wind region of a WR star is able to reproduce 
the spectrum (Waxman et al.|2007 l. However, the properties 
of the underlying WR star and its wind would be peculiar ( ,G] 
|2007l l. 

The absence of the jet break is consistent with the wide 
outflow produced without the pressure support from the sur- 
rounding star. Both the supernova shock breakout and the 
tidal disruption scenario produce wide outflows. 

Strong early optical/UV emission is consistent with com- 
ing from the jet front. The non-relativistic outer shells 
producing early optical/ UV blackbody em ission violate the 
source energetics ( |Ghisellini et al. [2007b i within the super- 
nova shock breakout model. In addition, the total radiated 
energy lO''^^ erg of GRB060218 is somewhat higher than ex- 
pected for the shock breakouts ( Colgate 1974 ) in general. 

The host of GRB060218 is a dwarf galaxy, and the dwarf 
galaxies are expected to have the highest tidal disruption rates. 
Only the dwarf galaxies and not the bigger galaxies are ex- 
pected to have the IMBHs in their nuclei, and thus disrupt the 
WDs. We crudely estimate the central BH mass 1 x IQ*Mq 
based on the stellar mass lO^^M© of the host. The supernova 
position is consistent with the center of the host dwarf galaxy, 
as it should be for the tidal disruption scenario. Supemovae 
unrelated to the central BHs generally have large displace- 
ments from the galactic centers. The supernova SN2006aj has 
the low estimated ejecta mass Mgj = 1 -2M0 consistent with 
the heavy WD. 

In sum, the tidal disruption model can explain all features of 
the source and only the scale of the jet power remains some- 
what arbitrary. Note that the BH mass 1 x 10* Mq is esti- 
mated in three entirely independent ways! A conclusion that 
GRB060218 was a tidal dismption by the IMBH in a nucleus 
of a dwarf galaxy at a distance 150 Mpc indicates that the 
IMBHs are abundant in the local Universe. 
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